Superdiffusion in the topological metal 
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We develop a non-perturbative theory to study large-scale quantum dynamics of Dirac particles 
in disordered scalar potentials (the so-called "topological metal"). For general disorder strength 
and carrier doping, we find that at large times, superdiffusion occurs. I.e., the mean squared 
displacement grows as ~ tint. In the static limit, our analytical theory shows that the conductance 
of a finite-size system obeys the scaling equation identical to that found in previous numerical 
studies. These results suggest that in the topological metal, there exist some transparent channels 
- where waves propagate "freely" - dominating long-time transport of the system. We discuss the 
ensuing consequence - the transverse superdiffusion in photonic materials - that might be within 
the current experimental reach. 
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The past years have witnessed that graphene [H [2] and 
topological insulators become an exciting frontier of 
condensed matter physics. A characteristic feature com- 
mon to these novel materials is the emergence of the 
massless Dirac particle. For example, it appears on the 
surface of three-dimensional strong topological insulators 
such as Bi2Se3 and Bi2Te3 [3H5]. This particle displays 
helical spin (that accounts for the sublattice structure in 
graphene) structure in momentum space and is "relativis- 
tic" : the energy dispersion is linear. On one hand, the 
discovery of the Dirac particle in these materials is fos- 
tering conceptual developments in physics, notably the 
Klein tunneling in condensed matter systems [6 ; on the 
other hand, with many striking properties, the Dirac par- 
ticle may find considerable potential applications in new 
electronic devices. 

For experiments and practical applications, the inter- 
play between Dirac particles and disorders is a key issue 
[7j. A complete theory must include the description of 
quantum dynamics of the Dirac particle in disordered en- 
vironments. This is a subject largely unexplored. Yet, 
there have been increasing evidences indicating that rich 
dynamic phenomena might occur in this system. Indeed, 
in the absence of disorders, experimental and theoreti- 
cal studies (see Refs. O [9] and references therein) have 
shown that the Dirac particle has already exhibited inter- 
esting dynamic behavior. In addition, being in the same 
symmetry class notwithstanding, a disordered Dirac par- 
ticle system and a normal metal with disordered spin- 
orbit coupling have transport properties of profound dif- 
ferences [TUHIl]. It is thereby conceivable that in the 
presence of disordered potentials, the interplay between 
wave interference and the relativistic energy dispersion 
may lead to even more interesting dynamic behavior and 
eventually, to unusual electric and optical properties. 

The purpose of this Letter aims at a systematical and 
analytical investigation of this subject. Specifically, we 
shall focus on the two-dimensional system with a single 
Dirac valley and subject to scalar disordered potentials 



F(x), the so-called "topological metal" pTO^. For this 
system, the quantum dynamics is characterized by a two- 
component spinor, VE*, obeying (We set h = 1.) 

idt^^H^, H=-iv{a''da: + aydy) + V{x.). (1) 

Here, v is the Fermi velocity and , i — 0, x,y, z stands 
for the Pauli matrices. The effective time-reversal sym- 
metry, i.e., ayH*ay — H, brings the system to the 
sympletic symmetry class. Below we develop a non- 
perturbative theory to study quantum dynamics ([T]) at 
large scales. 

Armed with this analytical theory, we find that for 
general disorder strength and carrier doping, the quan- 
tum transport of this system exhibits certain "anoma- 
lies". I.e., (i) at large times, quantum superdiffusion oc- 
curs: the mean squared displacement - characterizing the 
expansion of wave packets - is given by 



(5x2 (t) 



tint. 



(2) 



where p^^ is the density of states at the energy ep. It sug- 
gests that such a topological metal behaves as a "transi- 
tion" from Ohmic (Jx^ ~ t) to perfect (Jx^ ~ t^) metals. 
The t In t behavior, though of pure quantum origin (as we 
will see shortly) , is akin to a classical transport anomaly 
discovered long ago [T3]. Experiences in studies of clas- 
sical superdiffusion [T3] then suggest that the finding of 
quantum superdiffusion might have far-reaching impacts 
on many scientific branches including statistical physics, 
condensed matter physics, nonlinear dynamics, and pure 
mathematics. We notice that sufficiently away from the 
Dirac point, this expression is universal: it does not de- 
pend on disorders at all, as p^p converges to its clean 
limit. In contrast, near the Dirac point, this expression 
is non-universal because p^p is governed by disorders, (ii) 
In the static limit, our analytical theory shows that for a 
disordered sample of size L, the conductance g{L) follows 
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the scaling equation, 



irrespective of g. This is consistent with numerical find- 
ings [21 [T0HI2] . Remarkably, unlike a normal metal with 
disordered spin-orbit coupling, this system does not ex- 
hibit two-dimensional Anderson transition. The nature 
of anti-localization in topological metals is currently un- 
der intense investigations [10, 14J. 




FIG. 1: Multiple scattering of a Dirac particle by disordered 
potentials (pink). Inset: Klein tunneling. 

Let us first discuss a possible physical picture underly- 
ing Eqs. ^ and ^ and give their "intuitive" derivation. 
In doing so, we will see that (i) long-time transport of 
the system is dominated by transparent channels where 
waves propagate "freely" ; and (ii) the full quantum dif- 
fusion coefficient (including multiple scattering effects) 
displays logarithmic singularity. 

We shall focus on the Dirac point {ep = 0). In this 
case, the Klein tunneling [S] generally occurs (Fig. [l] in- 
set): upon hitting a potential barrier, a particle may be 
converted into a hole and subsequently undergoes per- 
fect transmission; similar scenario may happen to a hole. 
On general grounds, we expect (Fig. [T]) that assisted by 
the Klein tunneling, particles acquire a significant prob- 
ability of passing "freely" through disordered potentials 
before deflecting from the incidence direction, leaving a 
"free" flight path (much) longer than the distance be- 
tween two nearest disordered potentials. In the extreme 
case, they keep moving along the incidence direction and 
pass freely through all the disordered potentials encoun- 
tered, leaving an infinitely long free flight path ("trans- 
parent channels"). Thus, multiple scattering leads to a 
length distribution of the free flight path. 

The diffusion coefficient scales as ~ w^7^^, where the 
scattering rate 7 is inversely proportional to the length 
of the free flight path defined above and thereby includes 
multiple scattering effects. Averaging v'^j~^ with respect 



to the scattering rate distribution P(7) gives the long- 
time asymptotic of the diffusion coefficient, D{t~^ — )■ 
0) - d7P(7)7-i. Because of P(0) - (p.^w^^-i^ 
the integral over 7 suffers logarithmic divergence and is 
dominated by 7 ~ 0. That is, long-time transport of the 
system is dominated by the (almost) transparent chan- 
nels, i.e., 7 — > 0. Since the free flight path cannot be 
longer than vt at finite time, i.e., v^~^ < vt. the scat- 
tering rate has a lower bound, t"^ . (The upper bound 
is of minor importance as we are interested in the time- 
dependent behavior.) Taking this into account, 

D{t-^ ^ 0) - ^2^(0) / ^77-1 - Ini/Pep. (4) 

Substituting it into the mean squared displacement, 
(5x2 (t) _ ^(^-1 _^ o)t, we reproduce Eq. For 
a finite system, scales as ~ . In combination 
with the Einstein relation, Eq. ^ gives the conductance 
g{L) ^ e^lnL. Therefore, we reproduce Eq. ^ also. 
Surprisingly, Eq. ^ is insensitive to the details of ^(7). 
This universality may reflect the topological nature of 
the system. 

We are now ready to present the analytical theory 
and some technical details of the proof. We begin 
with weakly disordered scalar potentials that are Gaus- 
sian distributed, t/(x) = , V{x)V{yi') = ^'^(x - 
x'). Here, the dimensionless parameter go character- 
izes the disorder strength, and the overline stands for 
the disorder average. To proceed further we introduce 
the two-particle correlation function, (/)^(x — x',w) = 
^2il*^''^^^fF+;.^(^'^')G'^p(x',x), n = 0,x,y, where the 
matrix retarded (advanced) Green function obey {e—H± 
«(5)G'f'^(x,x') = (5(x - x')ct°. Notice that the disorder 
average restores the translational invariance and correla- 
tion functions thereby depend only on x— x'. As such, we 
may pass to the Fourier representation, (/>'^(x — x',aj) — > 

(/)'^(q,w) = 0^(g) and likewise, Gf'^(x,x') -> / d^{x- 

x')e-*P (''-''')Gf"*(x,x') = 5f'^(p). Here, Gf^"^ satis- 
fies (e - wa>a - Ef"4(p))5f '-^j-pj ^ ^0 ^ Pqj. ^g^^j^ 

disorders, one may invoke the self-consistent Bonn ap- 
proximation (SCBA) to calculate the self-energy, Sf'^, 
obtaining Sf^^ = 2^^"^^^^'^- (^^otice that the trace 
operation includes the momentum index.) The den- 
sity of states = — ifrlmC/^ is related to the imagi- 
nary part of the self-energy via = — '^ItrlmS^. Im- 
portantly, the SCBA gives an exponentially small en- 
ergy scale ~ e~'^'^ ^° [16]: for energies far below this 
scale (the Dirac regime), the density of states saturates, 

^ e""^^ while for energies far above it, converges 
to its clean limit ~ |e|. 

Then, by adaption of the method of Refs. [T71 [TH], it 
can be shown that in the hydrodynamic limit, g — >■ 0, 

a;0O(q)-q„r(g) = -p,,. (5) 
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This is the macroscopic continuity equation reflecting the 
particle conservation law, where 4'^{q) is the macroscopic 
density and {(j)^{q), (t>^{q)) the macroscopic current. The 
latter follows a Fick-like law, read 

r{q) = -iD{q)q^^\q) (6) 

in the Fourier representation. Here, D(ci) = v'^/(—iLu + 
j{q)) is full quantum diffusion coefficient, where 7(g) is 
the so-called "current relaxation kernel" given by 

liq) = - + l^^tT {q^a"Ag,U,Agyqpa^} . (7) 

Here, AC/,(p) ee tj,^+^(p + q) - Qt^ip), the elastic scat- 
tering rate ^ — — 2trlml]^ — v'^p^^/go, and f7g(p,p') is 
the two-particle irreducible vertex function. This equa- 
tion, together with the expression of D{q), builds up 
a bridge between the macroscopic hydrodynamics and 
the microscopic quantum dynamics. It is evident from 
Eq. ([7]) that localization effects - encoded by C/q(p,p') - 
introduce renormalization of the elastic scattering rate. 
In the absence of quantum interference, Eq. ([T]) leads to a 
Boltzmann diffusion coefficient D = v^t (with the classi- 
cal transport mean free time twice larger than the elastic 
scattering time). Eqs. ([5])-([7| constitute the framework 
of subsequent analysis. 

The closed set of macroscopic equations ([s]) and 
^ gives an exact relation between the optical con- 
ductivity, cr((?), and the full quantum diffusion coeffi- 
cient, D{q). To this end, let us introduce the den- 
sity response function defined as x(x — x',f — t') = 
i([*t(x, t)5'(x, t), ^-t (x', t')]) , with (•) the ground 

state average. Upon passing to the Fourier representa- 
tion, x(x — x', < — t') — x{q)i one may follow the deriva- 
tion of Ref. [T7] to show x{q) ^4>^{q) + P<LF- I'^ com- 
bination with the solution to Eqs. ^ and ([6|, it gives 

a{q) EE -e2 hm f ^x(<z)) = e^Pe^Diq)- (8) 

Thus, we justify the general Einstein relation for disor- 
dered Dirac particles. 

The system's behavior at large scales is determined 
by the current relaxation kernel, 7(9). Therefore, we 
proceed to perform non-perturbative analysis of this 
kernel dominated by the infrared divergences of the 
two-particle irreducible vertex function C/g(p,p'). Gen- 
erally, these divergences include two types: the singlet 
component of diffuson and cooperon |19j . It can be 
shown that the effective time-reversal symmetry leads to 
a mathematically rigorous theorem: 

The current relaxation kernel 7(g) suffers no diffuson- 
type infrared divergences. 



It generalizes the Vollhardt-Wolfle theorem discov- 
ered for spinless electron systems [17) . 

Thanks to this theorem, to calculate the two-particle 
irreducible vertex, Uq{p, p'), we need to consider the dia- 
grams composed of (the singlet component of) cooperon. 
For large frequencies, the dominant contribution to 
Uq{p,p') consists of single cooperon, giving a quan- 
tum correction to the (bare) elastic scattering rate ^ 
/ iS^i-'^ + Dk'r' order of 0{g^'). This is 
the well-known weak anti-localization, with the negative 
sign reflecting the destructive interference in the presence 
of the TT Berry phase [101 HH HSl |20] . It should be noted, 
however, that even for weak disorders, this quantum cor- 
rection is divergent in the low-frequency limit w — > 0. In 
this case the elastic scattering rate (and thus the Boltz- 
mann diffusion coefficient) undergoes strong renormal- 
ization and as such, the perturbative expansion in l/go 
breaks down. To go beyond the perturbation, we must 
take into account the full two-particle irreducible vertex 
for Eq. Q. As a matter of the reciprocity principle, 
Uqip,p') - <^0(p-l-p',c.) ^ [-tuj + DiLu){p + p')^]-\ 
where D{uj) = _D(q — 0, w). This amounts to the replace- 
ment of the Boltzmann diffusion coefficient in the weak 
anti-localization correction by the dynamical diffusion co- 
efficient D{uj). In doing so, we obtain a self-consistent 
equation of the dynamical diffusion coefficient, 

D 1_ fjFk 1 

that is expected to hold for arbitrarily low frequencies. 
Eq. ([9]) implies the validity of one-parameter scaling 
in topological metals, consistent with non-perturbative 
studies based on field theories |14j . 

Let us discuss below the physical consequences of the 
general theory described by Eqs. ([5])-(|9|. 

a) Superdiffusion. — We prepare a wave packet with the 
energy near ep {> 0) and let it evolve. Its expan- 
sion is characterized by the full time-dependence of the 
mean squared displacement given by (5x^(t) — / ^^(l — 
e~'"^*)D(w)/cLi^. For low frequencies, ujt <^ e^^^^3°, 
Eq. Q may be solved analytically, giving D{w)/D w 
^^2g^ In ■ Substituting this result into the expression 

of 5y?{t), we obtain Eq. ^ for t > re^'^'f". Noticing 
that ioj ^ t^^, we have D{t^^)/D w In |- namely 

Eq. Q. It is thereby justified that at large times, lo- 
calization effects lead to strong renormalization of the 
elastic scattering rate. 

b) Optical conductivity. — According to the Einstein re- 
lation for low frequencies, ujt ^ e~^'^^^°, the opti- 
cal conductivity reads <T(q = 0,lu) « (^)^ In that 
strikingly, is universal: it depends on the elastic scatter- 
ing rate logarithmically. 

c) Static conductance of a finite system. — We have fo- 
cused on the dynamic property of a bulk (infinite ex- 
tended) system so far. Now, we switch to the static con- 



ductance of a finite system of size L. This subject is 
currently under intense investigations [2j fT0Hl2] . In the 
static limit, Eq. (|9| is simplified to 



D 

m 



thanks to 



d^k 1 



7r.9o 7i<|k|<i (27r)2 



= 1 



(10) 



0. 



— r Kj. This equation gives the static 
diffusion coefficient £)(0) = D[l ' ^ 



^n— In— 1 



To translate I?(0) to a static size-dependent conduc- 
tance (identical to the conductivity in two dimensions), 
g{L), we apply a static, infinitesimally small electric field 
E to the system say, along the a:-direction. In response to 
the electric potential, {L—x)eE, a density profile is estab- 
lished across the sample, read ~eE dx'x{x~x' , 0)(i — 
x') — eEp^p{x — L). In obtaining this result, we have 
used the fact that in the real space, the static density 
response function is x(a; — x', a; = 0) = p^p5{x — x'). Re- 
sults from this density profile a uniform diffusive current, 
—e^p^pD{0)E, that is balanced by a macroscopic electric 
current, g(L)E. It then follows g{L) = pepD{0), giv- 
ing Eq. (|3| (where the conductance is in unit of /h.) 
For large samples, L/vt ^ e^^ ^° ^ the conductance 
g{L) = In that depends logarithmically on the 
elastic scattering rate and has a universal prefactor. 
d) General disordered potentials. — Qualitatively, the 
macroscopic equations ([5|, ([6| and ([s]) are the results 
of particle conservation, while Eq. Q reflects the nature 
of one-parameter scaling and one-loop self-consistency. 
As such, we expect them to be universal for a large class 
disordered potential. Speciflcally, the disorder potentials 
- remaining of the scalar type - may be strong, i.e., 
go ^ 1, and (or) exhibits long-ranged Gaussian corre- 
lations. These dramatic modifications notwithstanding, 
they merely affect the microscopic parameters namely 
Pep and D in Eqs. ([s]), ([6|, ([8| and Then, we may re- 
peat the discussions of a) — c), obtaining Eqs. ^ and (|3|. 
The latter is consistent with numerical findings pi fTHHT^ . 
despite that the method here differs essentially from the 
routine one based on the Kubo formula [TUj . 

In realistic electronic materials, electron-electron in- 
teraction interplays strongly with wave interference ef- 
fects (as observed in experiments on topological insula- 
tors [7j). Therefore, we discuss here the possibility of 
observing superdiffusion in photonic materials that are 
free of such an interplay (Fig. [2j left). We adopt the 
method that leads to the experimental discovery of the 
so-called transverse localization [3T] . In the present con- 
text, one may invoke the optical induction technique to 
fabricate a two-dimensional photonic crystal of thickness 
z that displays the honeycomb structure in the transverse 
(x-y) plane while is uniform in the longitudinal (z) direc- 
tion. As such, the photonic band structure exhibits two 
non- equivalent Dirac valleys, around K and K' (Fig. [2] 
right). One may further introduce random fluctuations of 





FIG. 2: Transverse superdiffusion. Left: a laser beam (blue) is 
launched into a two-dimensional disordered photonic crystal. 
The latter is uniform in the longitudinal (2) direction while in 
the transverse (x-y) plane, it displays a periodic background 
of honeycomb structure (green) with smooth random dielec- 
tric fluctuations (pink). At a sufficiently large distance z, the 
broadening of the probe beam scales as ~ zlnz. Right: the 
first Brillouin zone of the honeycomb lattice. 



refractive-index in the transverse plane that are smooth 
on a scale much larger than the lattice constant. Then, 
a monochromatic laser beam is launched into this disor- 
dered photonic crystal and undergoes diffraction broad- 
ening. At a distance z from the input plane, the intensity 
profile of the probe beam is measured. 

The probe beam excites only the Bloch modes in cer- 
tain Dirac valley, and has an initial width much larger 
than the scale over which the refractive-index fluctu- 
ates. Its propagation in the disordered photonic crys- 
tal is described by the (rescaled) paraxial Hemlholtz 
equation [3T], idzA = [— ^(9^ + dy) — n{x,y)]A, where 
A{x, y, z) is the envelop of the electric field smoothly 
varying in the longitudinal direction. Being static (hav- 
ing no-time dependence) notwithstanding, this equation 
possesses a perfect analogy to the Schrodinger equation, 
with z playing the role of "time" and the (negative of) to- 
tal refractive-index, ~n{x, y), of "potential". As a result, 
the diffraction broadening of the probe beam, A{x, y, z), 
mimics quantum wave function at time z. Because the 
"potential" is smooth enough, the inter-valley scatter- 
ing is suppressed. Effectively, this "quantum dynamics" 
is reduced to that described by Eq. Q, with the spin 
degree of freedom accounts for the sublattice structure. 
Thus, we expect that the broadening of the probe beam 
(in the transverse plane), JJ dxdy{x^ + y'^)\A{x,y, z)\'^ , 
scales as ~ zlnz for sufficiently large thickness z. We 
term this phenomenon transverse superdiffusion. 

In summary, we develop a non-perturbative theory of 
large-scale quantum dynamics of Dirac particles in dis- 
ordered environments. We show that for general disor- 
der strength and carrier doping, (i) a bulk topological 
metal exhibits superdiffusion ^ at large times; and (ii) 
the static conductance of a finite-size system follows the 



5 



scaling equation (|3|, consistent with the numerical re- 
sults [21 [T0HI2] . We discuss the ensuing consequence 
- the transverse superdiffusion in two-dimensional pho- 
tonic materials - that might be within the current exper- 
imental reach. 
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